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O ■ ABSTRACT 

ON 

o . 

\0 . Transitions of nuclear compositions in the crust of a neutron star induced by stellar 

G> ■ 

■ spin-down are evaluated at zero temperature. We construct a compressible liquid- 

Qh' drop model for the energy of nuclei immersed in a neutron gas, including pairing and 

O I shell correction terms, in reference to the known properties of the ground state of 

^ ', matter above neutron drip density, 4.3 x 10^^ g cm~^. Recent experimental values 

and extrapolations of nuclear masses are used for a description of matter at densities 
' below neutron drip. Changes in the pressure of matter in the crust due to the stellar 

^ I spin-down are calculated by taking into account the structure of the crust of a slowly 

and uniformly rotating relativistic neutron star. If the initial rotation period is ~ ms, 
these changes cause nuclei, initially being in the ground-state matter above a mass 
density of about 3 x 10^'^ g cm~^, to absorb neutrons in the equatorial region where 
the matter undergoes compression, and to emit them in the vicinity of the rotation 
axis where the matter undergoes decompression. Heat generation by these processes 
is found to have significant effects on the thermal evolution of old neutron stars with 
low magnetic fields; the surface emission predicted from this heating is compared with 
the ROSAT observations of X-ray emission from millisecond pulsars and is shown to 
be insufficient to explain the observed X-ray luminosities. 

Subject headings: dense matter nuclear reactions, nucleosynthesis, abundances 

stars: interiors stars: neutron stars: rotation X-rays: stars 
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1. INTRODUCTION 



A neutron star crust consists of a lattice of nuclei embedded in a roughly uniform 
neutralizing background of electrons and, at densities above neutron drip, in a sea 
of neutrons. So far, properties of the ground-state matter in the crust, such as the 
nuclei present and the equation of state, have been well investigated on the basis of 
laboratory data on nuclei and many-body calculations of the properties of uniform 
nuclear matter (see, e.g., Pethick & Ravenhall 1995). If the crustal matter, initially 
in its ground state, is compressed or decompressed, the matter is expected to deviate 
from nuclear equilibrium; the constituents of matter, giving its lowest energy state, 
now turn into another ones. It is expected, furthermore, that the matter being out of 
nuclear equilibrium undergoes nuclear processes such as electron captures, decays, 
neutron absorption and emission, and nuclear fusion and fission; these processes can 
release nuclear energy to lower the total energy of the matter. 

One of the astrophysical situations in which the pressure of matter in the crust 
is changing is in evolved neutron-star-black-holc binary systems; there may occur the 
tidal disruption of neutron stars by black holes. Lattimer et al. (1977) have investigated 
the compositions of expanding, initially cold, neutron star matter appropriate to these 
systems. It was noted that the compositions in the decompression bear a resemblance 
to those of hot r-process matter following a supernova explosion. Another situation is 
when a neutron star accretes matter from a companion star. If matter of ~ O.OIMq 
is accreted on to the surface of the star, the crust is composed totally of this accreted 
matter. Sato (1979) and Haensel & Zdunik (1990a) have examined the compositions of 
cold matter initially of a density of about 10^ g cm^"^ being compressed up to ~ 10^'^ g 
cm~^, which we may encounter in accreting neutron stars. Their conclusion was that 
the charge number of nuclei in the matter decreases as a result of electron captures 
and that, once neutrons drip out of the nuclei, the mass number also decreases by 
neutron emission. At the highest densities where Coulomb barrier between the nuclei 
is sufficiently low, pycnonuclear fusion was found to occur. 

ROSAT observations of X-ray emission from rotation-powered pulsars give us im- 
portant information on the surface emission from isolated (not accreting) neutron stars 
of various ages (Ogelman 1995). As the rotation is slowing down, the matter in the 
crust far away from the rotation axis is expected to be compressed due to decrease 
in the centrifugal force; on the other hand, the matter in the polar region may be 
decompressed due to increase in the length of the star along this axis. In this paper, 
we shall then investigate transitions of the nuclear compositions in the crust under 
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such compression and decompression and effects of the release of nuclear energy on the 
thermal evolution of neutron stars. 

During the time scale for the spin-down r^, which is > 10 yr as deduced from a 
magnetic dipole braking, neutrino cooling processes work so effectively that the tem- 
perature of the crust is below ~ 10^ K (see, e.g., Nomoto & Tsuruta 1987). For such 
temperatures, thermal energies are smaller than typical excitation energies of the sys- 
tem except in the outermost layers of the star where the thermal motion of electrons 
and ions affects the matter properties (Pethick & Ravenhall 1995). In describing matter 
in the crust of mass densities p > 10^ g cm~^ considered below, therefore, we assume its 
temperture to be zero. We also assume that this zero-temperature matter is in nuclear 
equilibrium before the compression and decompression proceed due to the spin-down. 

Nuclear processes can occur in the cold crustal matter during the spin-down only 
when two conditions are satisfied. The first is that the time scale for the processes 
must be less than r^, and the second is that the Gibbs free energy of the system should 
be lowered by the processes. We first consider these criteria for electron captures, f3- 
dccays, and neutron absorption and emission. The time scale for neutron processes 
occurring above neutron drip density via the interactions of neutrons in continuum 
states with nuclei (Lattimer et al. 1977) is negligibly small compared with Tg. Ac- 
cording to the estimates made by Bisnovatyi-Kogan & Seidov (1970) and Lattimer et 
al. (1977), the time scale for /9-processes, dependent on the nuclear structure and on 
the electron Fermi energies, may possibly be smaller than in the density region of 
interest to us, although this is still uncertain for the neutron-rich nuclei present near 
and above neutron drip density. The energy criterion for the occurrence of the pro- 
cesses considered may be re-expressed as the chemical potential of electrons (dripped 
neutrons) being across the energy threshold which is given by the difference in nu- 
clear masses between before and after /5-processes (neutron processes). Generally, this 
energy threshold depends strongly on pairing and shell effects in the nuclei, so that 
these effects should be duly taken into account in describing the nuclei in the crust. 
Next, let us inquire when nuclear fusion and fission can occur. The time scale for 
these processes is sensitive to the nuclear charge Z, since Z determines the size of 
Coulomb barrier between the nuclei and that of fission barrier (Lattimer et al. 1977; 
Sato 1979). Energetically, fusion and fission reactions can take place if the nuclear 
binding energy per nucleon becomes larger after these reactions. For both conditions, 
the occurrence of fusion and fission requires sufficiently small and large values of Z, 
respectively, compared with those for the equilibrium nuclei. 

For the purpose of examining the chemical evolution of the neutron star crust as 
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the star spins down, we construct a compressible liquid-drop model of nuclei in the 
inner crust (at densities above neutron drip) based on that of Baym, Bethe, & Pethick 
(1971; hereafter referred to as BBP) and use a description of matter in the outer crust 
(at densities below neutron drip) by Baym, Pethick, & Sutherland (1971; hereafter 
referred to as BPS). For the matter above neutron drip, not only is a nuclear surface 
energy term of BBP modified but also a shell correction term is added to the BBP 
expressions in such a way that the resulting compositions of the ground-state matter 
are consistent with those obtained by Negele & Vautherin (1973; hereafter referred to as 
NV) in their Hartree-Fock calculations of complete unit cells with an effective nucleon- 
nucleon interaction. A pairing correction term obtained by Moller & Nix (1992) as 
an empirical formula is also added to our nuclear model; an isospin dependence of 
nuclear pairing gaps found for laboratory nuclei was incorporated into their formula. 
Recent experimental values of the nuclear masses (Audi & Wapstra 1993) and their 
extrapolations (Moller et al. 1995) are built in the BPS expressions for the energy of 
matter below neutron drip. The resultant equation of state of the ground-state matter 
in the crust is found to be consistent with that of BPS. The next step is to calculate 
changes in the pressure of matter in the crust due to the stellar spin-down. We first 
describe the evolution of structure of the crust by using an equation of structure of 
a slowly and uniformly rotating relativistic star advanced by Hartle (1967) and the 
equation of state of BPS as well as by assuming a magnetic dipole spin evolution. 
Displacements of matter in the crust during the spin-down are then analyzed in the 
framework of general relativity. By comparing these displacements with the evolving 
pressure profile of the crust, we evaluate the variations in the pressure of the matter. 
Using the constructed nuclear model and the evaluated pressure variations as well as 
recalling the energy and time-scale conditions for the occurrence of nuclear processes 
described in the preceding paragraph, we investigate nuclear processes occurring in 
the crust during the spin-down, the resulting heat generation, and its effects on the 
thermal evolution of neutron stars. 

In § 2, we first construct a model for zero-temperature matter in the crust, and 
we thereby calculate the compositions and equation of state of the matter in nuclear 
equilibrium. In § 3, a model for the crust of a neutron star that is rotating down 
is described; changes in the pressure of the crustal matter induced by the spin-down 
are calculated therefrom. In § 4, we evaluate transitions of the compositions of the 
zero-temperature matter, undergoing the changes in its pressure as calculated in § 3, 
from the equilibrium ones obtained in § 2. In § 5, a heating rate associated with these 
transitions is estimated, and its relation to the ROSAT observations of millisecond 
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pulsars is examined. Conclusions are finally given in § 6. 



2. MODELS FOR THE COLD CRUSTAL MATTER 



In this section we construct an expression for the Gibbs free energy of zero-tempera- 
ture matter in the neutron star crust. We then investigate the equilibrium nuclei 
present and the equation of state for these nuclei by minimizing the Gibbs free energy 
of the system under charge neutrality. We assume that the matter is composed of 
a single species of nucleus with nucleon number A and proton number Z at a given 
pressure P. 



2.1. Inner Crust 



A compressible liquid-drop model originally developed by BBP is useful for de- 
scribing higher-energy states of cold matter in the inner crust on a footing equal to its 
ground state. BBP gave the energy Wn{A, Z, VN,Tin) of a spherical nucleus of volume 
V/v immersed in the uniform neutron gas of number density 

Wn = [(1 - x)mnc^ + xrripc^ + W(k, x)]A + Wc,ui(Z, Vn) + W^^f{A, Z, Vn, n„) . (1) 

Here mnimp) is the neutron (proton) rest mass, W{k,x) is the energy per nucleon for 
bulk nuclear matter of number density n = 2A;^/37r^ = A/Vn and of proton fraction 
X = Z/A [see equation (3.19) of BBP], FFcoui(-2^) ^n) is the Coulomb energy of protons 
distributed uniformly in the nucleus, 

Wcoul (2) 
5 rN 

with rN = {3VN/4:7iy^^ and W^^^{A, Z, Vn , Un) is the nuclear surface energy. 



^BBP ^ <7mkr„0)-W{k,x)] (^^_nn^ '/\2/S (3) 

with kn = {3n^nn/2y/^, a = 21.0 MeV, and uj^ = 16.5 MeV. 

In the ground state of matter above neutron drip density, we take note of the nuclear 
charge obtained by NV in their Hartree-Fock (HF) calculations of complete unit cells. 
The Z values calculated by BBP agree fairly well with those by NV at mass densities 
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p between 4.3 x 10^^ and 10^^ g cm^^, but significantly exceeds them with increasing 
density for p > 10^^ g cm~^ (see Figure 6 of NV). Such behavior of Z is due to 
overestimation of the nuclear surface energy by BBP, which was clarified by Ravenhall, 
Bennett, & Pethick (1972; hereafter referred to as RBP) through comparison with 
their HF calculations of the surface energy with a Skyrme interaction. We thus affix 
to H^surf^ ^ coefficient A acting to reduce this surface energy: 

f 1 J Tin = 

^=ltanh(g) , n„>0, 

where 

^^^^w{K,o) + lkn^^^^ (5) 

is the chemical potential of dripped neutrons (not including the rest mass), and Ci 
is the adjustable parameter to be determined below. Hereafter we replace W^surf^ 
equation (1) by IVsurf defined as Wgurf = ^W^t' and refer to the equation thereby 
modified as equation (1)'. 

The nuclear shell effects, giving rise to greater stability of the nuclei with neutron 
or proton magic numbers, also affect the charge number of the equilibrium nuclei; these 
effects were appropriately taken into account by NV through a proton spin-orbit force. 
We then add to equation (1)' a shell correction term that was originally developed for 
isolated nuclei by Myers & Swiatecki (1966) and extended to nuclei in the dripped- 
neutron regime by Sato (1979): 

■ F{N, N,)x + F{Z, Zi)x + F(Z, Z;)(l - x) 



W^shell = 5.8 



(^/2)2/3 

-0.26[(1 - x)x + I MeV (6) 



, M,_i < X < M, , (7) 



X^T—^-TTZY' (8) 

1 + exp 



with 
F(X,M,) = ^ 



where = A—Z is the number of neutrons in the nucleus, = {2, 8, 14, 28, 50, 82,126, 
184,---} is the set of neutron magic numbers for isolated nuclei, Zi = {2, 8,14,28,50,82,114,- 
•■} is the set of proton magic numbers for isolated nuclei, Z'- = {2, 8,20,28,40,50,82,114,- 
••} is the set of proton magic numbers for nuclei in the external neutron gas found by 
NV, and 

l^n = ^[^^(^' ^) + ^surf (^, Z, Vn, n„)]z,y^,„„ (9) 
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is roughly the chemical potential of neutrons in the nucleus. The factor x ensures that 
equation (6) reproduces the empirical shell formula of Myers & Swiatecki (1966) for 
isolated nuclei, and that the dripped neutrons smear the neutron shells and turn the 
proton magic numbers Z-i into Z[ (see Sato 1979). 

The parameter C2 has been determined together with Ci in such a way as to 
reproduce the Z values of the ground-state matter calculated by NV at eleven mass 
densities, so that Ci = 5.5 MeV and C2 = 0.1 MeV. Strictly speaking, we have 
obtained Z = 40 in place of Z = 32 (NV) at the highest density. In this determination 
we have evaluated the values of Z by minimizing the Gibbs free energy of the system 
under charge neutrality, as will be discussed in § 2.3. 

In Figure 1 we have compared the resultant surface tension. Eg = lVsurf/4vrr^, with 
the BBP and RBP results. The present values approach the HF values of RBP with 
increasing density (decreasing proton fraction) for p > 10^^ g cm~^ {x < 0.26) and 
reproduce the BBP values in the density range 4.3 x 10^^ g cm"'^ ^ P ^ 10^'^ g cm~^ 
(0.26 < X < 0.31). An improvement of the present surface energy over that of BBP 
has been thus clarified. 

We next introduce the nuclear pairing effects, stabilizing the nuclei with even neu- 
trons and even protons relative to the other types of nuclei, into the nuclear model 
hitherto determined by adding a pairing correction term obtained empirically by MoUer 
& Nix (1992) in their macroscopic model: 



W ■ 

^ ' pair 



(even-neutron, even-proton) 

r^^^-tyoilP- (odd-neutron, even-proton) 

r_^m^Q-iym&'^ (even-neutron, odd-proton) 

T^^^-t^^P ^ r^m^^-t^^S^ _ _|m|_ (odd-neutron, odd-proton) 



(10) 



with /iMN = 6.6 MeV, tmn = 5.55 MeV, and tym — 6.07. Here 5 = 1 — 2x is the relative 
neutron excess, and Bg is the surface area of a deformed nucleus divided by that of a 
spherical one. We simply set Bg — l. It is to be noted that they determined Wpair from 
experimental mass differences in the isospin range —0.1 < b < 0.25, and that we shall 
extrapolate VFpair into the dripped-neutron regime (0.4 < 5 < 0.7). We also note that 
this pairing correction term, being measured relative to the even-even nuclei, produces 
no change in the charge number of the equilibrium nuclei. 

We finally obtain the Gibbs free energy of the system above neutron drip per nu- 
cleon, as BBP did, according to 

9 = 1 (11) 

'T'Ar^cell 

with 

Exot = nNm^{A, Z, Vat, ^at, ej^n^ + (1 - nNV^)en{n^ , (12) 
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where Etot is the total energy per unit volume, which is in turn related to p as Etot = pc^, 
and njv is the number of nuclei per unit volume. Here, 



rriN = Wn{A, Z, Vn, n„) + Wl{Z, Vn, un) + W^sheii(^, Z, V^, n„) + W^^M, Z) (13) 

is the sum of the mass of a nucleus given by equations (1)', (6), and (10) and the lattice 
energy in the Wigner-Seitz approximation. 



with Tc = (3/47rniv)^/3^ 

Se = + + 1)^/^ - Ht + it' + 1)^/^} (15) 

with t — hiSn^UeY^^ / rrieC is the energy of the free electron gas per unit volume, 

£n = nn[W(kn, 0) + m^c^] (16) 

is the energy of the external neutron gas per unit volume, and 

Aceii = A + n„ f — - T/jv) (17) 

is the nucleon number in a spherical Wigner-Seitz cell, which is equal to the ratio of 
the number density of baryons rib to njv- 



2.2. Outer Crust 

To construct a model of matter in the outer crust, we follow a line of arguments of 
BPS and set the Gibbs free energy of the system per nucleon as 

(^«) 

with 

-Etot = nNmN{A, Z, hn) + e^irie) , (19) 
where is given by equation (15), and 

ruN = Mn{A, Z) + Wl{Z, njv) (20) 
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is the sum of the mass of a nucleus Mn{A, Z) and the energy of a bcc lattice, 

Wl = -0.895929 . (21) 

Here the nuclear mass, not having the atomic-electron binding energy subtracted out, 
is related to the atomic mass M^(A, Z) as M^(A, Z) — Ma{A, Z) — Zmf.c^. The values 
of Ma{A, Z) used here are taken from the experimental mass table of Audi & Wapstra 
(1993), when present. The masses of the remaining atoms are taken from the droplet 
(FRDM) mass table of MoUer et al. (1995). We note that the matter model described 
above is the same as that of Haensel & Pichon (1994) except some minor differences 
concerning the electronic contribution to binding energies, the electron exchange ener- 
gies, and the electron-screening effects. 



2.3. Ground-State Properties 



It is necessary to know the ground-state properties of matter in the crust for deter- 
mining the initial model for the equation of state and the nuclear compositions. The 
equilibrium conditions for the matter arc obtained by minimizing the Gibbs free energy 
per nucleon, given by equations (11) and (18) at fixed P under charge neutrality. 



Zn 



N ■ 



(22) 



For the matter in the inner crust, as investigated by BBP, the conditions for determining 
riN and Vn at fixed P, A, Z, and n„ are the pressure equilibrium condition. 



dm 



N 



dV, 



N 



p(G) ^1 . dW{kn-,Q) UnTlN OniN 



dkn 



1 - VnUn dun 



(23) 



A,Z,Vff,ni\f 



and the thermodynamic relation between P and rib, 



p = p(G) + n, 



^drie KrieJ ^ drtN 



(24) 



Z,Vn 



The remaining conditions for determining A, Z, and at fixed P, V^, and n^v are the 
optimization of A in a nucleus, 

d /P^Coul + W^surf + W^L- 



dA V 



x,nNA,njsiVN,nn 



<Sa 



(25) 
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with 

W^^M -hZ) + W^sheii {a - 1, Z, ^VW, ^n^r, n„) 



S 



A 



A-1 

Wp^riA Z) + W^sheii(^, 1/jv, njv, 



(26) 



the stabihty condition, 

Wz<lJie< Wz-i (27) 

with 

Wz = mjv(^, ^, Vat, un, Un) - mN{A, Z + 1, Vn, un, n^) , (28) 

where = dSe/dne is the electron chemical potential including the rest mass, and the 
equilibrium of the dripped neutrons with the neutrons in a nucleus, 

DA<f^n + rrinC^ < Da+1 (29) 

with 

DA = mN{A,Z,VN,nN,nn)-mN(A-l,Z,VN,nN,nn + - — ) • (30) 

V l — nNVNy 

In the conditions (27) and (29), the quantities Wz, Wz-i, Da, and Da+i correspond to 
the energy thresholds necessary for /3— decay, electron capture, neutron emission, and 
neutron absorption, respectively, to occur in a cell containing the nucleus {A, Z) . We 
note that inequalities (25), (27), and (29) are only necessary conditions for determining 
the compositions of the ground-state matter; hence, we should find the parameters A, 
Z, and Un, which give the minimum value of g, among those satisfying these inequalities. 

For the matter in the outer crust, a complete description of the minimization of g 
is given by BPS. In this case, the equihbrium conditions (24) and (27) remain useful 
if we set P^^^ = and eliminate the parameters Vn and n„. 

The neutron drip point for the ground-state matter has been calculated by minimiz- 
ing g given by equation (18) and finding the case in which g — rrinC^, as evaluated by 
BPS. The result indicates that a neutron begins to drip out of ^^^Kr at p = 4.34 x 10^^ 
g cm-^ P = 4.93 X 10-^ MeV fm^^ and //g = 26.21 MeV. The agreement with the 
results of BPS and of Haensel & Pichon (1994) has been thus confirmed. 

Using the equilibrium conditions discussed above, we have obtained the equilibrium 
nuclear compositions for eleven cases of densities below neutron drip as well as for eleven 
cases of densities above neutron drip, which were chosen by NV, and determined the 
equation of state for these compositions; the result has been hsted in Table 1. We have 
confirmed that the calculated nuclear compositions (including the chemical potentials 
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He and fin) are consistent with those of Haensel & Pichon (1994) in the outer crust as 
well as those of NV in the inner crust. The equilibrium nuclear compositions listed 
in Table 1 will be used as the initial model for the nuclear compositions in § 4. As 
depicted in Figure 2, the equation of state in the present calculations agrees graphically 
with the equation of state of BPS which includes that of BBP in the dripped-neutron 
regime. We shall thus use the equation of state of BPS as that consistent with the 
present result for estimating the structure of the crust in § 3. 



3. SPIN-DOWN AND CHANGES IN PRESSURE OF THE CRUSTAL 

MATTER 



In this section we first construct a model for the crust of a zero-temperature neutron 
star rotating at a given angular velocity. By using a magnetic dipole approximation 
for the rotational evolution, the evolution of structure of the crust is determined. We 
finally estimate the displacements of elements of matter in the crust during the spin- 
down and then the corresponding changes in their pressures. We use units in which 
c = G = 1 in § 3. 



3.1. Model for the Crust of a Rotating Neutron Star 



As a first step towards the estimates of the changing pressures of matter in the 
crust during the spin-down, we consider a pressure profile of the crust of a neutron 
star that is in hydrostatic equilibrium and is rotating at a uniform angular velocity 
Q = 27r/Prot) where Pj-ot is the rotation period as seen by a distant observer. In a 
non-rotating equilibrium configuration, we assume the stellar mass M and radius Rs 
to be 1.4 Mq and 10 km. We fix the total number of baryons in the star for arbitrary 
Q. General relativistic effects on the crust model may be deduced from the following 
quantity included in the redshift factor at the surface of the non-rotating star: 

2GM 



RsC 



2 



0.4 . (31) 



A fully general relativistic treatment of the configuration of the crust is thus required. 
If the star is rotating with Prot = 1 ms, being similar to the shortest period 1.6 ms 
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among those observed from pulsars, the ratio between the typical magnitude of the 
rotational energy of the star and that of its gravitational potential is roughly 



2GM ^''^ 
We can then consider the rotation as a small perturbation on a non-rotating configu- 
ration. Consequently, we find that general-relativistic equations of structure of slowly 
and rigidly rotating stars obtained by Hartle (1967) exactly up to order are useful 
for our purpose of determining the pressure profile of the crust of a rotating neutron 
star. 

Such a pressure profile may be obtained by determining the transformation of a 
surface of constant pressure in the equilibrium configuration of the crust of a neutron 
star rotating at a given O from a spherical one of radius R in the non-rotating equilib- 
rium configuration subject to conservation of the number of baryons in the star. We 
denote the position of the surface of constant pressure in the rotating configuration 
by the ordinary polar coordinates r and 6 in which the polar axis is taken to be the 
rotation axis. At a given 9, the radius of the surface may be written as (Hartle 1967) 

r^R + ^{R,e) + 0{n^), (33) 

where 

CiR,e)^CoiR)+UR)P2icos0) (34) 

is the displacement of order from its non- rotating position. Here Co{R) denotes the 
change in the mean radius of the surface, and ^2(-R)-P2(cos 9) represents the quadrupolar 
deformation of the surface, where P2(cos^) is the second-order Legendre polynomial. 
For the present purpose, we should then evaluate the radius it! of the surface of constant 
pressure in the configuration of the non-rotating neutron star crust and its displacement 
^(i?, 9) due to the rotation. 

We first calculate the pressure profile, P{R), of the crust of a non-rotating neutron 
star by using the equation of state of BPS in the Tolman-Oppenheimer-Volkoff (TOV) 
equations (see Shapiro & Teukolsky 1983): 



dR i?2 



(35) 



where 

e^(«) = 1 - ^ (36) 
R 



with 



m{R) = M - / ' (ir47rrV(r) (37) 

JR 
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is the redshift factor. The pressure profile P{R) obtained by solving equation (35) 
inward from the surface of the star is shown in Figure 3, together with the density 
profile ni){R). Here we have assumed that the crust dissolves into uniform nuclear 
matter at a density of ~ 0.1 fm~^, as suggested by Lorenz, Ravenhall, & Pethick 
(1993). 

We proceed to estimate the displacement, ^{R, 6), of the surface of constant pressure 
in the crust from the non-rotating one due to the rotation, using the results for this 
displacement obtained by Hartle (1967) and Hartle & Thorn (1968). Since the thickness 
of the crust is small compared with Rg as depicted in Figure 3, it can be assumed that 
^{R, 9) ^{Rs, 9) throughout the crust. To calculate the displacement of the surface 
of the star ^{Rs,9), we note the solutions of Einstein's field equations for a perfect 
fluid in the gravitational fleld of the rotating conflguration that is characterized by the 
metric having the following form up to order fl'^ (Hartle 1967): 



ds^ = -e'^Wfl + 2[ho{r) + /i2(r)P2(cose)]}(dx°)' 



+e 



^ ^ 2[mo(r)+m2(r)P2(cosg)] | 
r — 2m (r) 



+r^l + 2[v2{r) - /i2(r)]P2(cos 0)}{de^ + sin^ ^[^0 - uj{r)dx^]^} , (38) 



where z/(r) and m(r) are given by equations (36) and (37). At the surface of the star, 
the solutions for the first- and second-order quantities contained in equation (38) are 
obtained by Hartle (1967) as 



ho{Rs) = 



uj{Rs) = 



SM 



Rs (l 



h2{Rs) = 



J2 



2M \ 
Rs J 

5 



2J 



+ 



J2 



Ri (1 - ^) 



1 —] 

MRl V 8M3 



mo (P.) = 5M - — 



« - f) (§ 



m^iRs) = {Rs - 2M) 



6J^ 
Ri 



V2{Rs) = —fTA + 



MQ - J2 



-Q2 ('§ - 1 



(39) 
(40) 

(41) 
(42) 
(43) 
(44) 



R's 4M3^^^(^^_2M)"' VM 

where J is the total angular momentum of the star related to the moment of inertia 
/ of the star as J — IQ,, 5M \s the change in mass of the star from its non- rotating 
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value, Q is the quadrupole moment of the star, and Q2"' is the associated Legendre 
polynomial of the second kind. At fixed baryon number, 5M is given by (Hartle 1970) 

5M = -7Q2 . (45) 

The quantity of ^2(-Rs) can be derived from the integral of the equation of hydrostatic 
equilibrium for the rotating configuration [see equations (28), (89), and (91) of Hartle 
(1967)] as 

«.. = -{i^S^.M..)}|(i-f). (.e) 

The parameters which still remain to be determined are 7, Q, and ^o(-Rs)- These 
values, dependent on the properties of neutron star matter, have been taken on the 
basis of the results for the rotating configuration at a stellar mass of about 1.4 Mq 
obtained by Hartle & Thorn (1968), using the equation of state that is hard rather 
than soft, as 

I = 10^^ g cm-3 (47) 
Q = 0.04 (^)'Mi?2 (48) 

URs) = 0.06 i^—J Rs , (49) 

where is the angular velocity at Prot = 1 ms. Here we have taken into account 
differences between the rotations at fixed central density and the rotations at fixed 
baryon number, the former being assumed by Hartle & Thorn (1968). We have thus 
obtained 

= (-0.039 + 0.149 sin^^) f-^f . (50) 

Equation (50), as well as the solutions of the TOV equations obtained by using the 
equation of state of BPS, gives the present model for the crust of a neutron star rotating 
at a given Q. 

In Figure 4 wc have drawn the surface of the star derived from equation (50) 
ai Q = Qq together with the surface of radius Rg at Q = 0. Wc observe in this 
comparison that the rotation acts to shorten (enlarge) the surface of the star in the 
polar (equatorial) region. 
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3.2. Rotational and Structural Evolution 



We assume that the evolution of the rotation of the star is determined by the usual 
magnetic-dipole model for pulsars (see Shapiro & Teukolsky 1983) as 

^ "ft = 0) . (51) 

Here B is the constant magnetic-dipole field at the magnetic pole of the star which 
is taken to be on the equator, and t denotes the age of the star. The constant B 
is then expressed in terms of the rotation period Pj-ot and its derivative P^t as P ~ 
^/lQ^^P^^^P^t•, the values of B that are determined by the quantities Prot and Ppot 
observed from pulsars range from ~ 10^ G to ~ 10^^ G (Taylor, Manchester, & Lyne 
1993). By using equation (51), the time scale for the spin-down may be estimated as 

r, ~ 6 X 10^[10^^G/B{G)f[Qo/^Kt = 0)? Y^- 

Let us now ask how the structure of the crust is varying as the rotation of the 
star is slowing down according to equation (51). We assume that the total number 
of baryons in the star remains constant during the spin-down, and that the star is in 
hydostatic equilibrium and is rotating rigidly at t = 0. Since the period derivative Prot 
observed from pulsars is extremely small (typically 10~^^ s s~^), the configuration of 
the star can be considered to maintain its hydrostatic equilibrium and rigid rotation. 
Then, the displacement of the surface of the star given by equation (50) as well as the 
perturbative quantities given by equations (39)-(44) may be redetermined at a given t 
as ^{Rs,9,t), u{Rs,t), etc., by substituting fl(t) into Q included in each equation. In 
addition, we neglect the variation in the equation of state of matter in the crust during 
the spin-down, an assumption to be justified by calculating the equation of state of the 
matter deviating from its ground state in § 4. The evolution of structure of the crust 
is thus described in terms of the quantity ^{Rs,6,t) and of the pressure profile of the 
non-rotating configuration obtained by using the equation of state of BPS. 



3.3. Changes in the Pressure of the Crustal Matter 



For our main purpose of investigating transitions of the compositions of the cold 
crustal matter during the spin-down, it is essential to calculate the Lagrangian changes 
in pressure associated with elements of matter in the crust. Such Lagrangian changes 
in pressure may be obtained by combining the displacements of matter elements from 
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their initial positions with the Eulerian changes in pressure obtained in § 3.2. For 
convenience, we denote the initial position of a matter element by the depth Ar, 
the distance measured from the surface of the star along the radial line, and by the 
polar angle 6. We then express the position of this element at a given t by the depth 
Ar+r7(Ar, 9^ t) and the polar angle ^^ + a(Ar, 9^ t), as shown in Figure 5. The quantities 
7] and a correspond to the displacements in the radial and angular directions between 
the element and an observer who is initially located at the same point as the element 
and is comoving with the surface of constant pressure in the radial direction. 

The direction in which elements of matter in the crust move as the star is spin- 
ning down gives the relation between the quantities rj and a. Let us assume that the 
direction of movement of an element present at (r, 9, t) is parallel to the sum of the 
gravitational and centrifugal forces acting on the element. In the hydrostatic equi- 
librium configuration considered here, these forces are in balance with the pressure 
gradient. As can be seen from the radius of the surface of constant pressure given by 
equation (33), the 9 component of the pressure gradient is of order whereas the r 
component is of zeroth order. We thus find a to be of order Q"^, by comparison with 
1] of order Q^. Consequently, the quantity a can be neglected within the framework of 
the present perturbative treatment. 

To find the quantity ri{Ar, 9, t), we note that each element of matter in the crust is 
to be identified by counting the number of baryons inward from the surface of the star 
along the radial line. This identification may be made by following the baryon number 
conservation law (see Misner & Sharp 1964) according to 



dRsm9{rh-^[l + A{r,9,t^0)]] nb{R) 
= dRsm9{r^e-'^[l + A{r,9,t)]} nb{R) , (52) 

jRs-/\r-r]{^r,9,t) I J r=R+^{Rs,9,t) 

where nb{R) is the density profile of the non-rotating configuration shown in Figure 3, 

and ^ 

A(r,g,t)^ ^^^;'^'^ -1 (53) 
r^'e 2 sin 9 

is the factor of order determined by the determinant g of the metric of equation 
(38) and by the time component of the 4-velocity of the matter, 

u\r,9,t) = Q-'^ |l + ^rV('')[n(i) -c^(r,i)]^sin2^-/io(r,i) -/i2(r,i)P2(cos^)| 
+0(0^) . (54) 

We then set A{r,9,t) A{Rs,9,t) in equation (52), by noting the correspondence 
with the approximation, ^{R,9,t) ~ ^{Rs,9,t), available for thin crusts as discussed 
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in § 3.1. The values of A{Rs,9,t) can be calculated from equations (36), (38), (39), 
(41)-(44), (54) as 

A{Rs,e,t) = h-''^'''^Rl{[m-^{Rs,t)]''-u;{Rs,ty}sm''e 



+2[v2{Rs:t) - h2{Rs:t)]P2{cOs9) + 



mo{Rs,t) + m2{Rs,t)P2{cos6) 



Rs-2M 



-0.043 + 0.094 sin^^) 



(55) 



We have thus obtained the quantity rj{Ar, 6, t) by retaining only the terms of order fl^ 
in equation (52) as 



r){Ar,e,t) = 



e 2 



R=Rs-Ar 



R^UbiR) 

{f^s ^ ^(R) 

\A(R„ 0,t^O)- A(R„ 0, t)] / dRR^e'—ndR) 
jRs-Ar 

+ [^{R„9,t^O)-aRs,0,t)] dR 

JRs-Ar dR 



d 



R'e-— 



nb{R) 



(56) 



We now turn to the evaluations of changes in the pressure of matter in the crust 
induced by the spin-down. Recall the definition of ri{Ar,6,t) shown in Figure 5 and 
the evolving pressure profile of the crust described in § 3.2. Then, the change until an 
age t in the pressure of an element that is initially located at the point denoted by Ar 
and 9 may be written as 



6P{Ar,e,t) = P{R)\r=r,-a 
= -v{Ar,0,t) 



r-'q{Ar,e,t) 

dP{R) 



— P{R)\R=Rs-Ar 



dR 



+ O(Q^) , 



(57) 



R=Rs-Ar 



where P{R) is the pressure profile of the non-rotating configuration shown in Figure 
3. We have calculated the values of r]{Ar,9,t) according to equation (56), and, by 
substituting these values into equation (57) , we have obtained the values of 5P{Ar, 9,t). 
The result has been expressed as 



5P{Ar,9,t) , ,^^^^ = 0)^-^(^)2 



= (-0.11 + 0.35 sin^^)- 



(58) 



P{Rs - At) 

Here we have omitted the Ar dependence of 6P{Ar,9,t)/ P{Rs — Ar), since it makes 
negligible differences. 
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Figure 6 illustrates as a function of 9 the values of 6P{Ar,6,t oo)/P{Rs — Ar) 
given by equation (58) in which we set Prot{t = 0) = 1 ms; this quantity denotes 
the change in the pressure of an element of the crustal matter at an age of far later 
than Tg divided by its initial pressure. We thus see that the matter is compressed in 
the equatorial region of the crust and is decompressed in the polar region due to the 
spin-down, as noted in § 1. 



4. TRANSITIONS OF NUCLEAR COMPOSITIONS IN THE CRUST 



Let us consider nuclear processes induced in the cold crustal matter, initially in 
nuclear equilibrium, by the changes in its pressure due to the spin-down as discussed in 
the preceding section. We first ask how the matter being compressed or decompressed 
departs from its ground state in which the system satisfies the equilibrium conditions 
(23)-(25), (27), and (29). We note that the time scale for pressure readjustment, 
corresponding to the sound travel time across a nucleus (Lattimer et al. 1977), is 
negligibly small compared with the spin-down time scale r^, so that the system should 
always satisfy the pressure conditions (23) and (24). Due to the assumption that 
the system contains a single species of nucleus at a given pressure, we may replace 
the condition (25) by the condition that the number of baryons, Aceii, in a single cell 
remains constant. If nuclear fusion (fission) were to occur, however, this process would 
double (halve) the value of A^eii- The conditions (27) and (29) ensure stabihty of the 
matter against (3— and neutron processes, respectively. While one of these processes is 
in progress, the system, not satisfying the corresponding stability condition, is out of 
quasi-equilibrium. Here we note that the time scale for neutron processes is far shorter 
than Tg-i and we assume that the time scale for /3— processes is also smaller than r^. 
Then, the system may be considered to retain quasi-equilibrium during the spin-down. 
We shall thus pay attention only to the energy criterion for the occurrence of nuclear 
processes in analyzing the chemical evolution of the crust. 

By noting the Gibbs free energy of the system per nucleon given by equations (11) 
and (18) as well as the conditions to be satisfied in the compression and decompression 
as discussed above, we have examined transitions of the composition of a unit cell 
from the equilibrium one listed in Table 1 due to the change in its pressure during the 
spin-down as follows: (i) At each pressure shifting from the equilibrium value Peq, we 
have calculated the Gibbs free energy per nucleon under the charge neutrality (22) and 
the pressure conditions (23) and (24) for the initial nuchde {A, Z, A^eii) as well as for 



18 



the nuclides {A,Z±1, A,,n), {A ± 1, Z, A,,n), {2A, 2Z, 2Aeeii), and {A/2, Z/2, Aceii/2), 
resulting from /3-processes, neutron processes, fusion, and fission, respectively, (ii) By 
finding the case in which the Gibbs free energy for the nuclide arising from the reactions 
becomes equal to the energy for the initial one, we have determined the type of reaction 
expected to occur first and its threshold pressure Pth- (hi) At the pressure equal to Pth, 
we have determined the transitions to the next nuclides according to the criterion that 
the Gibbs free energy calculated for the nuclide present after the reactions under the 
conditions (22)-(24) should be lower than that for the previous one, until the system 
attains quasi-equilibrium. (iv) By following a line similar to (i)-(iii), we have examined 
a series of processes expected to occur with a further change in pressure. In (i)-(iv) 
we have assumed that, if f3- and neutron processes were energetically possible at the 
same time, neutron processes would go first. 

In addition to the changes in the nuclear compositions, we have evaluated the nu- 
clear energy q thereby released per nucleon according to the thermodynamics at zero 
temperature as g = Qi — Qf-, where g.i is the Gibbs free energy per nucleon for the 
nuclide before the reactions, and gj is that for the next nuclide. As far as f5— and 
neutron processes are concerned, not only the reactions which release no energy but 
those which release some can take place. If the Fermi energies of electrons (dripped 
neutrons) become equal to the energy thresholds for (3— (neutron) processes during 
the compression or decompression, the electrons (neutrons) at the sea levels disappear 
or come out without supplying the medium with energy. Since these reactions proceed 
in a quasi-equilibrium condition, we shall refer to such processes as quasi-equilibrium 
processes. On the other hand, the reactions can liberate some energy when the Fermi 
energies of electrons (dripped neutrons) are in excess of the energy thresholds for elec- 
tron captures (neutron absorption), i.e., > Wz-i {l^n + irinC^ > Da+i), as well 
as when those are below the energy thresholds for /3— decays (neutron emission), i.e., 
He < Wz {fin + iTT'nC^ < Da)- In the former case, the electrons (dripped neutrons) with 
energy between Wz-i {Da+i) and /Xg (/U„ + ninC^) are captured by the nuclei, which 
leads to excitation of the nuclei and to creation of holes in the electron (neutron) 
Fermi sphere. In the latter case, the nuclei add the electrons (neutrons) with energy 
between ji^ {j^n + nT-nC^) and Wz {Da) to the electron (neutron) Fermi sea and simulta- 
neously undergo excitation. In both cases, therefore, the reactions reduce the system 
to a non-equilibrium state; hereafter we shall denote such processes by non-equilibrium 
processes. At zero temperature considered here, the following relaxation towards quasi- 
equilibrium results in liberation of the energy qAcdi per unit cell; for neutron absorp- 
tion, neutron emission, electron captures, and /3— decays, qAceii ^ fJ^n + ^nC^ — Da+i, 
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qAccW ^ Da- jJ^n- rrinC^, qAccW ^ jJ^e- Wz-i, and qAcdi ^ Wz - /J^e, respectively. The 
relaxation processes involved, the resulting heat generation, and the applicability of 
the zero-temperature approximation will be discussed in § 5. 

We proceed to show the results for the transitions from the equilibrium nuclides 
to the other ones during the compression of matter in the equatorial region of the 
crust due to the spin-down. The increase in its pressure can amount to about 24% for 
Prot{t = 0) = 1 ms as shown in Figure 6. In Figure 7 we have drawn the diagrams 
that show the changes in nuclides induced by such or a little larger increase in pressure 
from six equilibrium nuclides present at mass densities pcq above neutron drip. 

In the initial density region 4.3 x 10^^ g cm~^ < pcq ^ 3 x 10^^ g cm~^, as shown 
in Figures 7a-7c, electron captures occur first when P = Pth, at which the electron 
chemical potential fie, increasing due to the compression, attains the energy threshold 
Wz^i- These processes are quasi-equilibrium ones. Thereafter, a succession of neutron 
emission and electron captures takes place at the same pressure as non-equilibrium 
processes. This result indicates that the nuclei undergoing non-equilibrium neutron 
emission (electron captures) have energy thresholds Da {Wz-i) larger (smaller) than 
fin + nin(? (/ie) in contrast with the initial nuclei. This contrast is partly because the 
initial nuclei have even proton and neutron numbers and hence receive greater stability 
from the pairing effects, which are described by equation (10), and partly because they 
have closed shells of protons leading to local energy minima, which are determined by 
equation (6). A quasi-equilibrium state, in which the system includes the even-even 
nuclei in most cases, follows this series of reactions and lasts until the pressure arrives 
at the threshold for the forthcoming electron captures. 

For initial densities of about 3 x 10^^ g cm~^, neutron absorption begins to take 
place along with /3-processes as a result of the increase in the chemical potential of 
dripped neutrons during the compression, as can be seen from Figure 7d. This result 
implies that, at such initial densities, the nuclei are immersed in the neutron gas dense 
enough for the threshold pressure for neutron absorption to become comparable to 
that for electron captures. For higher initial densities up to about 2 x 10^^ g cm^^, as 
shown in Figures 7e and 7f, neutron absorption rather than electron captures occurs 
first as quasi-equilibrium processes when P = Pth at which /i„ + m„c^ ^ Da+i- These 
processes are accompanied by neiitron absorption which occurs at the same pressure 
in a non-equilibrium condition and hence releases some nuclear energy. The resulting 
state in which the system contains the even-even nuclei remains in quasi-equilibrium, 
until further increase in pressure leads to the next two-stage neutron absorption. The 
two-step feature of the reactions stems from the fact that the energy thresholds Da+i 
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depend on the parity of the reacting nuclei due to the nuclear pairing effects; this 
feature can be observed schematically from Figure 8a. 

At initial densities below neutron drip, it is shown that electron captures take place 
first at the threshold pressure higher than Pgq by more than 40%. Such increment in 
pressure, however, is larger than that expected from the compression due to the spin- 
down. This result suggests that the energy thresholds for electron captures, Wz-i, 
formed mainly by the closed shells of neutrons (protons for Ni) and pairing gaps in the 
equilibrium nuclei present in this density regime (see Table 1) are large enough to keep 
the nuclei stable against the reactions during the spin-down. 

As can be observed from the results mentioned above, no fusion reaction is allowed 
energetically by the compression during the spin-down in the whole density region of 
the crust. This is because electron captures, leading to a reduction in the nuclear 
charge, do not proceed sufficiently for the occurrence of fusion. 

We now summarize the nuclear processes expected to occur in the compressed 
matter during the spin-down and the energies thereby released. By comparing the 
threshold pressure Pth for the processes occurring first (see Figure 7) with the pressure 
to be increased by the spin-down with Prot(^ = 0) = 1 ms in the equatorial region 
of the crust (see Figure 6), we find that, when Prot(^ = 0) ^ 1 iiis, almost all the 
processes occurring in the compression are neutron absorption by the nuclei initially 
present in the density region 3 x 10^^ g cm^^ < pcq ^ 2 x 10^"^ g cm^'^. In Table 2 
we have shown the changes in the nuclear compositions in such a density regime with 
increasing pressure up to 24% and the energy q released per nucleon. We see that both 
the ratio of Pth to Pcq shown in Figure 7 and the energy released by non-equilibrium 
neutron absorption as listed in Table 2 tend to decrease as the initial density (proton 
fraction) increases (decreases), a feature coming primarily from the isospin dependence 
of Wpair given by equation (10). 

For the purpose of examining the change in the equation of state of the compressed 
matter from the equilibrium one during the spin-down, we have also calculated the 
equation of state for the compositions of matter which are altered from the equilibrium 
ones listed in Table 1 by an increase of 30% in pressure. Figure 9 depicts the resulting 
equation of state, together with the equilibrium equation of state listed in Table 1 and 
with the equation of state of BPS. We observe in this figure that the non- equilibrium 
equation of state, as well as the equilibrium one, agrees graphically with that of BPS in 
the nearly whole density region. Such agreement between the non-equilibrium and BPS 
equations of state may justify our neglecting effects of the variation in the equation 
of state during the spin-down on the structure of the crust, as noted in § 3.2. It is 
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interesting to notice the deviation of the non-equiUbrium equation of state from the 
equiUbrium one found for densities between 4.3 x 10^^ and 3 x 10^^ g cm~^. Such 
deviation originates from the fact that the threshold effects prevent neutron emission, 
leading to softening of the equation of state in this density range, from proceeding 
enough for the number density of dripped neutrons n„ to reach its equilibrium value 
at the same pressure, as suggested by Haensel & Zdunik (1990b). 

We conclude this section by showing the results for the changing compositions of the 
decompressed matter in the polar region of the crust during the spin-down. We have 
found that, basically, the processes occurring in the decompression are inverse processes 
against those caused by the compression in the whole density regime. This result is due 
to the fact that the electron and external neutron gases, controlling the occurrence of 
13- and neutron processes, undergo a reduction in their chemical potentials during the 
decompression in contrast with the case during the compression. For Prot{t = 0) = 1 
ms. Figure 6 shows that the pressure of the matter in the vicinity of the rotation axis 
decreases by about 11% during the spin-down. In consequence of such decrease in 
pressure, neutron emission is found to take place as two-stage reactions composed of a 
quasi-equilibrium process and a non-equilibrium one at initial densities between 5x 10^'^ 
and 2 x lO^"' g cm^"^. Figure 8b depicts a schematic of nuclear levels relevant to these 
two-stage reactions. In Table 3 we have shown the changes in the nuclear compositions 
in such a density regime with decreasing pressure down to 11% and the energy per 
nucleon released by neutron emission. In comparison with Table 2, we observe that the 
changes in nuclides during the decompression and those during the compression from 
the same equilibrium nuclides are analogous with each other in terms of the energy 
release and of the pressure intervals between the adjacent thresholds for the reactions. 



5. HEAT GENERATION AND NEUTRON STAR THERMAL 

EVOLUTION 



The energies released by neutron absorption and emission that occur during the 
compression and decompression turn completely into heat through relaxation processes 
such as scattering between neutrons in the continuum states and 7 de-excitation of the 
reacting nuclei. The heat thereby deposited locally in the substance is ~ 0.1 MeV 
per reaction, which is in turn carried away by electrons through their scattering with 
phonons, impurities, and other electrons (Flowers & Itoh 1976). The time scale for the 
heat conduction across a unit cell can be estimated as Tff ~ cprl/n, where cp is the 
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specific heat per unit volume at constant pressure, which is simply assumed to be the 
sum of that for the degenerate nonrelativistic neutrons of density n„ and that for the 
degenerate relativistic electrons of density rig, and k, is the thermal conductivity. We 
then note that the time interval, tj, between when the pressure reaches the threshold 
for the reaction in a cell and when in the neighboring cell is ~ TsTc/ i], where rj is given 
by equation (56). At typical densities and temperatures of the substance around the 
reacting cell, cp ~ 10^° erg cm~^, n ~ 10^*^ erg cm^^ s~^, Vc ~ 10 fm, and 
?7 ~ 10 m, so that th <^ t/. We can thus assume that the heat production by the 
reaction in a cell makes no change in the energy criterion for the reactions occurring 
next in the surrounding cells, which has been designed at zero temperature in § 4, and 
hence in the results shown in Tables 2 and 3 for the chemical evolution of the crust. 

The total heat, Qtot, deposited in the star during the spin-down with Prot{t = 0) = 1 
ms may be estimated from the following: (i) The energy per unit volume released by a 
non-equilibrium neutron process that occurs in the matter initially present at densities 
between 3 x 10^'^ and 2 x lO^'^ g cm^^ (0.02 fm~^ ^ ni, < 0.1 fm"^) is evaluated from 
Tables 2 and 3 as njg ^ 4 eV fm~^, irrespective of densities, (ii) The reacting material is 
initially located at the interior of depths Ar between 450 and 700 m, as can be seen from 
Figure 3. (iii) The number of the non-equilibrium processes, N{Ar), occurring during 
the spin-down in a unit cell that is initially located at the depth Ar on the equatorial 
plane is estimated from Table 2 and Figure 3 as N{Ar) ^ 15(Ar — ri)/(r2 — ri), 
Ti < Ar < r2, where ri = 450 m and r2 = 700 m. (iv) The ratio between the number 
of the non-equilibrium processes taking place on the cone of fixed 6 and that on the 
equatorial plane is ~ sin^|l — 1.46cos^^| as deduced from equation (58). By using 
(i)-(iv) and by taking account of the redshift effect on energy generation inside the 
non-rotating star (Thorn 1977), we have roughly estimated Qtot as 



This value of Qtot is found to be smaller than the initial rotational energy of the star, 
/Qo/2, by six orders of magnitude; therefore, Qtot has neghgible effects on the rotational 
evolution determined by the magnetic dipole braking. 

We now consider a heating rate H{t), the heat deposited in the star by the non- 
equilibrium neutron processes that occur per unit time at an age t. Note that both 
the pressure intervals between the adjacent thresholds for the non-equilibrium reactions 
occurring in a unit cell and the energies thereby released take on almost constant values, 
as can be seen from Tables 2 and 3. Then, we may assume that the heat generated 




10**^ erg . 



(59) 
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until a given age is proportional to the variation in the pressure of the reacting material 
due to the spin-down until that age as described by equation (58). We thus obtain 



W -f-"Wr (60) 

"0 



Consequently, the heating rate may be expressed as 
H{t) fa 5 X 10-'^E{t) 
5 X 10^2 



S(G) 



1010 G 



2 n(t=o)4 



1.54 



B(G) 



IQiOG 



^ erg s"^ , (61) 



n(f=0)2 f (yr) 1 ' 
lOSyr J 



where 

£;(i) = -m{t)n(t) (62) 

is the spin-down power derived from equation (51). 

To ascertain the practical utility of equation (61), we recall that the occurrence of 
neutron processes in the crustal matter requires the changes of at least a few percent 
in its pressure, as seen from Tables 2 and 3. In consequence, the condition for use of 
equation (61) is roughly estimated from equation (58) as 

Q(i = 0)2 - Q(i)2 > O.ml . (63) 

We thus see that the present heating mechanism works when Prot(^ = 0) < 3 ms and 

t > O.lTs. 

Figure 10 illustrates as a function of the age the heating rate given by equation 
(61) in which Prot(t = 0) = 1 ms and P = 10^ 10^ 10^°, 10i\ lO^^ q xhe age at 
which H(t) starts to decline remarkably corresponds to the spin-down time scale r^. 
We first consider the age of up to ~ 10® yr when neutrino emission from the interior 
of the star dominates over photon emission from the surface (see, e.g., Nomoto & 
Tsuruta 1987). In the light of models for the neutron star thermal evolution including 
some heating sources (see, e.g., Cheng et al. 1992; Reisenegger 1995; Van Riper, Link, 
& Epstein 1995), we expect that, at such times, the present heating source inside 
the stars even with high magnetic fields ~ 10^^ G does not contribute to the surface 
emission considerably compared with the contribution of the initial cooling of the stars. 
For t > 10® yr, on the other hand, it is generally predicted that the heat deposited 
initially is gone, and that the surface emission, dominant over the neutrino emission, 
is controlled entirely by some heat generation inside the star. It is noteworthy in 
Figure 10 that the heating rate for P ~ 10^ G lasts during a term far longer than 
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10^ yr. We thus see that the old neutron stars with low magnetic fields, characteristic 
of millisecond pulsars, are the most significant examples in investigating effects of the 
present heating mechanism on the neutron star thermal evolution. 

In Figure 11 we have depicted the bolometric and X-ray (0.1-2.4 keV) luminosities, 
L^i and L^, predicted from the heating rate (61) as a function of the spin-down 
power E, together with the X-ray luminosities obtained from the i^O^/ir observations 
of some millisecond pulsars in the energy range 0.1 — 2.4 keV (Banner, Kulkarni, & 
Thorsett 1994) and with an empirical formula for X-ray luminosities coming from 
the stellar magnetosphere as proposed by Ogelman (1995). Here we have evaluated 
the bolometric luminosity by assuming a complete balance between the photon 
emission from the surface of the star and the heat production inside the star as well as 
by taking into account the redshift effect on the surface emission in the configuration 
of the non-rotating star (Thorn 1977) as 

The X-ray luminosity has been then obtained by assuming a spectrum of the emis- 
sion from a spherical blackbody of radius Rs having the uniform surface temperature 
as observed by a distant spectator (Nomoto & Tsuruta 1987): 

where ctsb is the Stefan-Boltzmann constant. We have selected the pulsars satisfying 
the condition (63) among those chosen by Danner et al. (1994). 

We observe in Figure 11 that the X-ray luminosities of millisecond pulsars are far 
larger than L^. In comparison with the empirical result for magnetospheric X-ray 
luminosities obtained by Ogelman (1995), however, the X-ray luminosities observed 
are predicted to contain a large amount of nonthermal magnetospheric component in 
addition to a thermal component. We may thus conclude that the heat deposited in 
the star by neutron absorption and emission during the spin-down contributes little to 
the X-ray luminosities from millisecond pulsars; such a contribution is buried in the 
other components such as the magnetospheric one. 



6. CONCLUSIONS 

We have calculated the changes in the pressure of zero-temperature matter in the 
crust of an isolated rotating neutron star during its spin-down and the resulting transi- 
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tions of the compositions of the matter initially in nuclear equilibrium. We have found 
that the matter is compressed and decompressed in the equatorial and polar regions 
of the crust, respectively, as the rotation of the star is slowing down. The correspond- 
ing change in the equation of state has been shown to have negligible effects on the 
structure of the crust. Through the compression (decompression), the nuclei initially 
present at densities between 3 x 10^^ and 2 x 10^^ g cm~^ have been found to absorb 
(emit) neutrons, when Prot(^ = 0) < 3 ms. The heating rate arising from these pro- 
cesses has been estimated; we have clarified its marked effects on the thermal evolution 
of old neutron stars with low magnetic fields of ~ 10^ G. It has been shown that the 
X-ray (0.1-2.4 keV) luminosity predicted from the present heating mechanism is far 
smaller than the values observed from millisecond pulsars. 

The uncertainty affecting the occurrence of nuclear processes, which has been found 
in the dripped-neutron regime, stems primarily from the present extrapolation of the 
empirical formula for nuclear pairing gaps to higher relative neutron excess d, as noted 
in § 2, and possibly from the omission of pairing correlations in the external neutron 
gas. Fortunately, the pairing effects in the neutron gas, leading to stabilization of 
the continuum states of neutrons, would be insensitive to the energy thresholds for 
reactions, to which, however, the pairing correlations between nucleons inside a nucleus 
are crucial. We expect that the validity of the extrapolation used here will be confirmed 
by microscopic calculations of the pairing gaps in the nuclei far from the /3-stability 
valley with a realistic nucleon-nucleon interaction. 

Another problem which remains to be investigated is the detailed description of 
the transition from roughly spherical nuclei in the external neutron gas to uniform 
nuclear matter. Lorenz et al. (1993) and Oyamatsu (1993) discovered that the phases 
containing nuclei of nonspherical (cylindrical or planar) shapes may exist at densities 
lower than the point at which the matter becomes uniform. If so, these nuclear shapes, 
affecting the discrete states of nucleons inside a nucleus, may have nonnegligible con- 
sequences for the occurrence of nuclear processes and hence for the heat generation. 
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TABLE 3 



Nuclear Compositions in the Inner Crust during the 
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^ Initial equilibrium nuclides. 
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FIGURE CAPTIONS 



FIG. 1. The surface tension as a function of x, the proton fraction in nuclei. 

The sohd curve is the result obtained in the present calculations, the dotted curve is 
the BBP result, and the dashed curve is the RBP result. 

FIG. 2. Equation of state of the crustal matter in its ground state. The solid 

circles denote the present result, and the line is the BPS result. 

FIG. 3. Pressure and density profiles of the crust of a non-rotating neutron star 

of mass 1.4 Mq and radius 10 km, as calculated from the equation of state of BPS. 
The solid line is the pressure P (in MeV fm~^), and the dashed line is the number 
density of baryons nj, (in fm~^). The distances R and Ar (in km) are measured from 
the center and surface of the star, respectively. We have set the density at which the 
phase with nuclei changes into the liquid phase as 0.1 fm~^; the depth Ar at this 
phase boundary is represented by the dotted line Ar = r2, where r2 = 700 m. The 
dotted line Ar = ri, where ri = 450 m, corresponds to the depth at a density rib ~ 0.02 
fni-3. 

FIG. 4. Surfaces of a neutron star in the rotating and non-rotating configurations 

at fixed total baryon number. The sohd curve is the surface of the star rotating at a 
rotation period of 1 ms, and the dashed curve is the surface of the non-rotating star of 
radius 10 km. The z axis is taken to be the rotation axis. 

FIG. 5. Definition of the quantities r) and a, denoting the position of an element 

of matter in the crust at an age t (open circle). Its initial position (solid circle) is 
denoted by Ar and 9. The surface (a) is the surface of a neutron star rotating at an 
angular velocity fl{t). 

FIG. 6. The change in the pressure of an element of matter in the crust at an 

age when the rotation decays completely, being divided by its initial pressure, as a 
function of the polar angle denoting its position. The initial rotation period Prot(^ = 0) 
is taken to be 1 ms. 

FIG. 7. Changes in nuclides due to increase in pressure at initial densities above 

neutron drip. The solid circles arc the initial equilibrium nuclides at a density pcq and 
a pressure Peq- The crosses are the nuclides in the matter being in quasi-equilibrium. 
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The pressures Pth are the threshold ones for the reactions occurring first. 

FIG. 8. Schematic diagrams of nuclear levels relevant to two-stage neutron (a) 

absorption and (b) emission, which take place in the matter initially present above a 
density of about 3 x 10^^ g cm~^ during the compression and decompression, respec- 
tively, due to the spin-down with Prot(^ = 0) = 1 ms. When the Fermi energy of dripped 
neutrons including the rest mass, /x„ + m„c^, crosses the energy threshold, Da+i {Da), 
for neutron absorption (emission) by the even-even nucleus {A, Z) , the nucleus changes 
into the next one {A + 1, Z) [{A — 1, Z)] in a quasi-equilibrium condition. Successively, 
another neutron absorption (emission) occurs, which in turn releases the energy qAcew 
in the medium. The resulting state in which the system contains the even-even nucleus 
stays in quasi-equilibrium. The two-step feature of the reactions described above is 
attributable to the differences between the values of Da for odd A and those for even 
A, which occur due to the nuclear pairing effects. 

FIG. 9. Changes in the equation of state of zero-temperature matter in the 

crust from the ground-state one due to the compression. The solid circles show the 
equilibrium equation of state, and the crosses denote the equation of state of the matter 
deviated from nuclear equilibrium by an increase of 30% in its pressure. The solid line 
is the equation of state of BPS. 

FIG. 10. Heating rates as a function of the age of neutron stars with various 

values of surface dipole magnetic fields. The initial rotation period Prot(^ — 0) is taken 
to be 1 ms. 

FIG. 11. Photon luminosity from the surface of a neutron star as a function of 

the spin-down power. The solid line is the bolometric luminosity, L^^i, predicted from 
the heating rate H{t), and the dashed line is the X-ray (0.1-2.4 kcV) luminosity, Lx, 
predicted from H{t) by assuming a blackbody spectrum. For comparison, the X-ray 
luminosities observed from some millisecond pulsars (Banner et al. 1994) are shown. 
Triangles denote upper limits; asterisks detections; a diamond a possible detection. 
The empirical curve (Ogclman 1995; dotted line) for X-ray luminosities coming from 
the stellar magnetosphere is also plotted. 



33 



I I I I 




/ 

/ 

J I I I I I I I I I I I I I I I I I I I 



0.1 0.2 0.3 0.4 0.5 

x=Z/A 



FIGURE 1 
IIDA & SATO 



10 




109 lO'O IQll 10'2 1013 IQ14 



p (g cm-3) 



FIGURE 2 
IIDA & SATO 




FIGURE 3 
IIDA & SATO 




FIGURE 4 
IIDA & SATO 




FIGURE 5 
IIDA & SATO 




FIGURE 6 
IIDA & SATO 



45 



ISI 40 



35 



T r 



— I 1 1 1 1 1 1 

p^^=5. 012x10-4 MeV fm-^ 
-Peq=4. 667x1011 g cm-3 



1 



t(a) Pth/Pe,= l-41 
I I I I I I I I I 



125 130 
A 



45 



— I — I — I — I — I — I — I — I — I — 
P_= 1.483x10-3 MeV fm-^ 

eq 



35 



- p^q=2. 663x1012 g cm-3 
ISI 40 h 



(b) 



Pth/p..=i-32: 



Ih/ eq 

J I I I I I I I L 



130 



135 
A 



140 



55 



ISI 50 - 



45 



T 



T 1 1 r 



T 



T r 



P^q= 1.371x10-2 MeV fm-^ 
p^^= 1.497x1013 g cm-3 



1 



(c) 

_^ \ L 



Pth/Peq=l-23j 

J I I I I I 



190 



195 



55 



ISI 50 - 



45 



— I — I — I — I — I — I — I — I — I — 
P^^=4. 746x10-2 MeV fm-^ 
p^^=3.432xl0i3 g cm-3 



H n 



Pth/Peq=l-14j 
J I I I I I I I I 



(d) 



215 



220 
A 



225 



45 



I I I I I I I I I I I I I I I I I I I I I I I 
P^^=0.1777 MeV fm-^ 
-p^^=8. 010x10*3 g cm-3 



35 



(e) 



Pth/Peq=l-06j 



45 



I I I I I I I I I I I I I I I I I I I I I I I I I 
p^^=0.4094 MeV fm-^ 
- p^q= 1.334x1014 g cm-3 



N 40 — • — — ^ — ^ — ^ — — — ^ — N 40 — • — ^ — ^ — — — ^ — ^ — ^ — — — ^ — ^ 



35 



(f) 



p../p.„=i-03: 



th' eq 

I I I I I I I I I I I I I I I I I I I I I I I I I 



215 220 225 230 235 

A 



280 285 290 295 300 305 

A 



FIGURE 7 
IIDA & SATO 



(a) 





~ ^cell ^ 








Da+1 





D 



A+2 



(A, Z) 

even-neutron 
even-proton 



(A+1, Z) 

odd-neutron 
even-proton 



(A+2, Z) 

even-neutron 
even-proton 




(A,Z) — ► (A-1,Z) — ► (A-2,Z) 

even-neutron odd-neutron even-neutron 

even-proton even-proton even-proton 



FIGURE 8 
IIDA & SATO 




109 lO'O IQll 10'2 1013 IQ14 



p (g cm-3) 



FIGURE 9 
IIDA & SATO 



39 
38 
37 
36 
35 



00 34 
A 33 
K 32 

o 31 
^ 30 

29 

28 

27 



-B=10i2G 


P,,,(t = 0)=l ms_ 


- 10iig\ 


— 


- IQioG \ 




- lO^G 




- lO^G 







12 3 



4 5 6 7 
log t (yr) 



8 9 10 



FIGURE 10 
IIDA & SATO 




FIGURE 11 
IIDA & SATO 



